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^ ■ ABSTRACT 

We present a simple model based on the dark halo approach which provides a useful 
y—{ ■ way to understand key points determining the shape of the non-Gaussian tails of the 

' dark matter one-point probability distribution function(PDF). In particular, using 

C****" , the scale-free models with power-law profile of dark halos, we derive a simple analytic 

expression for the one-point PDF. It is found that the shape of the PDF changes at 
the characteristic value of 6* which is defined by the smoothed density of a halo with 
the characteristic mass M 4 at the epoch. In cold dark matter models with top-hat 
1 smoothing filters, the characteristic smoothed density at present time typically takes 

| the value <5* » 1 for a small smoothing scale R t h ~ IMpc/ft. and conversely 6* -C 1 

-»»« . for a large smoothing scale R t h > 10Mpc//i. On the range 6/6* < 1, the shape of the 

i-^ ' PDF is almost solely determined by the outer slope of halos and scales as a power- law. 

The resultant non-Gaussian tails of PDF then resemble the log-normal PDFs in that 
q ■ range and show a good agreement with N-body simulations, which can be ascribed to 

' the universality of the outer slope of the halo profile. In contrast, tails of one-point 

c/3 , PDF in the range 6/6* > 1 basically follow the steep exponential tails of the halo 

mass function, which exhibit a strong sensitivity to both the outer slope of the halo 
profile and the initial power spectrum. Based on these results, the discussion on the 
• <-h , PDF of galaxy distribution and the application to the weak lensing statistics are also 

/\i • presented. 

h ; 
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1 INTRODUCTION 

The probability distribution function (PDF) of the cosmic matter density is one of most fundamental quantities characterizing 
statistical properties of large-scale structures. In most standard scenarios of the structure formation, the large-scale structures 
are thought to form by gravitational instabilities from small Gaussian density fluctuations. A deviation in the PDF from 
the Gaussian distribution is thus most likely caused by the nonlinear gravitational evolution and is a key to understand the 
structure formation history. The higher order moments of the density field such like the skewness and kurtosis have been 
used to study the non-Gaussian nature in cosmic density field (e.g., Peebles 1980; Juszkiewicz, Bouchet & Colombi 1993; 
Bernardeau 1994a). The one-point PDF contains more statistical information than them. 

The galaxy PDFs have been measured from various catalogs (e.g., Hamilton 1985; Bouchet et al. 1993; Kofman et al. 
1994) and shows significant non-Gaussianity. Indeed it was pointed out that the galaxy PDF is well fitted by the log-normal 
distribution (e.g., Coles & Jones 1991; Kofman et al. 1994). The dark matter PDFs measured from TV-body simulations are 
also fitted by the log- normal distribution very well (e.g., Bernardeau & Kofman 1995; Taylor & Watts 2000). Although the 
galaxy distribution cannot be regarded as a direct tracer of the dark matter distribution because of unknown biasing relation 
between them, the similarity in their PDF suggests that there exists a nature common to their distributions. 

The origin of the phenomenological log- normal PDF has not yet been clarified. Bernardeau (1994b) and Benardeau & 
Kofman (1995) argued that the perturbative construction of one-point PDF shows a strong dependence on the primordial 
power spectrum and the PDF approaches the log-normal form only when the initial spectrum is proportional to P(k) oc k n 
with n = — 1. On the other hand, using a series of high-resolution N-body simulations, Kayo, Taruya & Suto (2001) found that 
the one-point density PDFs with Gaussian initial conditions are approximately described by the log-normal distribution in 
the strongly non-linear regime. Remarkably, the accuracy of log-normal model turns out to become almost insensitive to the 
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underlying power spectrum of density fluctuation, in contrast to the prediction from perturbation theory. While the results 
of Kayo et al. (2001) show a clear dependence on the primordial initial spectrum in a weakly non-linear regime, consistent 
with the perturbation theory, there remains crucial to clarify the reason why the PDFs in the simulations tend to become 
insensitive to the initial power spectrum and approach to the log-normal form in a strongly non-linear regime. 

The primary purpose of this paper is to discuss this issue using a simple analytic model. That is, employing the so-called 
dark halo approach as advocated recently, we construct a model for the non-Gaussian tails of the one-point PDF, which 
includes the essential ingredient for gravitational clustering and is, at least, consistent with N-body simulations. The basic 
assumption of our model is that the high density regions after smoothed over with a filter function are mainly attributed to 
a single massive halo. Clearly, this assumption is violated at the low density regions and we thus restrict our analysis to the 
highly non-Gaussian tails of mass density PDF. Then, the dark halo approach allows us to derive a simple analytic expression 
for the one-point PDF in cases with the scale-free initial spectra and power-law profiles of dark halo, from which we explain 
the qualitative behavior of tails of PDF and its dependences on the halo density profile and on the initial power spectrum. 
Further, adopting the realistic halo profile and the mass function, the results of our analytic model are quantitatively compared 
with N-body simulations with various initial power spectra. 

In section ^, the basic ingredients for the dark halo approach and an analytic model of tails of one-point PDF are 
described. Section]^ presents the detailed investigation of the dependence on the halo profile and the sensitivity to the initial 
power spectra on the non-Gaussian tails of one-point PDF. Assuming the scale-free initial spectrum and the power-law profile 
of dark halos, we analytically derive an approximate expression for non-Gaussian tails of PDF. Usingthis result, we attempt 
to clarify the sensitivity of the one-point PDF to the initial spectra and the halo profiles. In section W the analytic model of 
one-point PDF is quantitatively compared with those obtained from cosmological N-body simulations. Finally, section [| is 
devoted to the conclusions and the discussion. 



2 MASS DENSITY PDF FROM THE DARK HALO APPROACH 



2.1 Basic ingredients 

The basic idea of the dark halo approach which has become very popular recently (e.g., Seljak 2000; Ma & Fry 2000a; 
Scoccimarro et al. 2001) is that (i) all the dark matter in the universe belong to one of the virialized clumps (dark halos), 
and (ii) the matter distribution in the universe is approximated by a sum of the matter distribution within each halo. To 
construct an analytic model of mass density PDF from the dark halo approach, we need to specify the halo mass function 
and the density profiles of halos. In the following, we describe these two ingredients. 

For the halo mass function, the approximate prescription for counting the number of dark halos has been first presented 
by Press & Schechter (1974) and recently an accurate analytic fitting model has been proposed by Sheth & Tormen (1998). 
Both models are simply expressed in a unified manner as follows: 
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where p is the mean mass density and <5 C rit is the critical over-density of the spherical collapse model(see Nakamura & Suto 
1997 for useful fitting formula). The quantity a M denotes the linear rms fluctuations of the matter density in a sphere of radius 
Rm = (3M/47rp) 1//3 . The numerical coefficients a and p are set to a = 1, p = for the Press-Schechter mass function and 
a = 0.707, p — 0.3 for the mass function by Sheth & Tormen. The normalization factor A is determined by J n(M)MdM = p. 
Throughout the paper, we extensively use the mass function by Sheth & Tormen (1998). 

The density profile of dark halos is defined by a spherical average for each halo and is given by a function of radius and 
halo mass. The numerical simulation by Navarro et al. (1996,1997; hereafter NFW) revealed that the outer slope of profile 
becomes —3, while the inner part asymptotically approaches —1. In this paper, we adopt their fitting form: 

PNFw(r; M) = Sc 

P (r/r.)(l + r/r a ) 2 ' 1 ' 

where 5 C and r„ respectively denote the characteristic over-density and the scale radius. These parameters are not independent 
and are characterized in terms of concentration parameter c = r v i r /r s , with r v i r being virial radius of the halo. Using this 
quantity, the characteristic over-density 5 C is given by S c = (A v i r /3)c /[ln(l + c) — c/(l + c)], where A v i r denotes the virial 
over-density of spherical collapse model (Nakamura & Suto 1997). For the relation between the concentration parameter and 
the halo mass, we adopt the following fitting form according to the simulation by Bullock et al. (2001), 



c(M) 
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(3) 



where A/* corresponds to the characteristic mass scale satisfying a" M (Af*) = 5 C rit • Below in comparing the PDF with the 
N-body simulation (SecQ), the parameters c„ and (3 are chosen as (c»,/3) = (10,0.15) for cold dark matter(CDM) model 
and (10, (n + 3)/6) for scale- free model with spectral index n. The systematic influences of the uncertainty of c(M) on the 
evaluation of PDF are discussed in detail in Appendix B. 

Indeed, recent high-resolution N-body study reveals that the inner slope of halo profile might not converge to — 1, but 
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it rather approaches —1.5 (Moore et al. 1999; Fukushige & Makino 2001a, b; see also Jing & Suto 2000). Hence, it becomes 
important to clarify the profile dependence on the non-Gaussian tails of PDF. In this paper, in addition to the NFW profile, 
we further treat the halo models characterized by a single power-law profile: 

/Opowor(r;M) 

P 

which enables us to derive the analytic expressions of PDFs (see Sects) . Here, the characteristic density and scale radius are 



(~) > (0<«<3) (4) 

to derive the analytic expressions of PDFs (see 
determined by the relation between virial mass and virial radius 
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2.2 Mass density PDFs 

Provided the basic ingredients, we now describe a simple analytic model of mass density PDF. 

Let us first describe the cumulative distribution function P(> S;R), which gives a probability that a given density field 
smoothed over the radius R is larger than 8. Defining the smoothed halo profile: 

Pamooth (r;R,M) = / dV W{\r - r'\;R) p(r';M), (6) 



with the function W(x; R) being the spherically symmetric smoothing kernel. Since we only consider simple spherical symmetric 
density distributions of halos, one can calculate the characteristic radius r(S; R, M) by solving the following equation: 

p^ ooth [r(S;R,M);R,M] =5 + 1 ^ 5 (7) 
P 

Now we introduce the volume factor V(> 8; R, M) which represents, for one halo with mass M, the comoving volume in which 
the smoothed mass density field is always greater than 8. Since the smoothed halo profile is a monotonic decreasing function 
of the radius, the volume factor is simply given by, 

V(>S;R,M) = ^-{r(8;R,M)} 3 . (8) 

The volume factor (^) explicitly depends on the halo mass and thereby the volume fraction of the entire region greater 
than S is evaluated from the integration of equation (^) over the different halo mass weighting the comoving number density 
of dark halos. Indeed, this volume fraction is equivalent to the cumulative distribution function P(> <5; R). We thus obtain: 

P(> 8;R) = / dM V(> 8; R, M) n(M). (9) 



Once provided the cumulative distribution, it is easy to calculate the one-point PDF from the definition P(> 8;R) 
ft 



f+°° dSP{S;R): 
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Equation ( |Io| ) with equation (j^) is the heart of our subsequent analysis. It provides a simple analytic estimate for the 
degree of non-Gaussian tails, which is, in principle applicable to both weakly non-linear and strongly non-linear regimes. 
Note that the expression (^) heavily relies on the assumption that any of the dark matter particles reside in each dark halo. 
This would be true for high density parts of the mass density PDF, but might not be the case for the low density regions, 
— 1 < 8 ,6 1. We shall therefore focus on the high density tails, 8 ^> 1. As shown in the following sections, this approach is 
very useful to understand the essential physical mechanisms that determine the behavior of the high density tail of the PDF. 

While the above model prescription is general and does not restrict the choice of the smoothing filter, hereafter, we 
specifically treat the density field smoothed over the top-hat filter function: 

W th (r; fluO = —^3- 9(i? th - M), (11) 
where the quantity Q(x) denotes Heaviside step function. 



3 DEPENDENCES OF POWER SPECTRUM INDEX AND HALO PROFILE 

Before addressing the comparison with N-body simulations, it is instructive to consider how the halo density profile and 
initial power spectrum alter the non-Gaussian tails of PDF in a simple analytical manner. In this section, simply assuming 
the power-law model of halo density profile (^) and scale-free initial spectra P(k) oc k n , we here derive an approximate 
expression for the one-point PDF P(8;R t h)- 

To compute the PDFs, one must first evaluate the volume factor (0), which is given in terms of the characteristic radius 
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r(6; R t h, M). In Appendix A, the exact expressions for top-hat smoothed density profile is presented in the case of power-law 
profiles (eqs.]26|[^7|). The resultant smoothed profiles become constant at an inner region r <C Rth and show the power-law 
behavior at r 3> -Rth (see left panel of Fig.^J) . Based on these exact results, in Appendix A, we further obtain the approximate 
expression for the smoothed density profiles, which enable us to calculate the quantity r(S; R t h, M) analytically (see eq. |2q] ) : 



ooth(r;R th ,M) 



5o(Rt h ;M) ; r < (3=2) ,a R th 

(12) 

where the function So(Rth, M) represents the maximum value of the smoothed density profile, p smoot h(0; -Rth, M)/~p, and is 
given by a function of halo mass: 

/ r> \ —a y 7W- >. —a/3 

<5o(Rth,Af) = A vir (j^J =A vir (^J , (13) 

The quantity Mr means the effective halo mass inside the smoothing radius R t h, Mr = (47r/3)pA v i r Rf h . For a given 
smoothing length R t h, equation ([u^]) relates the mass of a halo with its characteristic smoothed density. Adopting the above 
approximation (|l^), one can easily solve (Q), which yields 

f ; 5>5 {R th ;M) 

r(5;R th ,M)={ . (14) 

[ j^MM} 1 ^^ . S <5o(R th ;M) 

Then, with a help of ([l3|), substitution of equation ( |l4] ) into (^|) leads to 

f ; 5>8 (R th ;M) 

V(>S;R th ,M) = { . (15) 

l ('*1T ^rf° ; 5<S (R th ;M) 

Clearly, the volume factor decreases as <5~ 3 / a and is cut off at the maximum density So{Rth', M). This indicates that for a 
certain value of S, the dominant contribution to the integral in cumulative distribution P(> 6; Rth) is mainly attributed to 
the halos whose maximum density is larger than that. 

Now, turn to focus on the evaluation of the integral in (^|). For the cosmological models with initial scale- free spectra, 
the rms fluctuation a M in (^) can be expressed as 

M \ -("+3)/6 

(16) 



rit (m») 



Substituting the mass function (|lj) and the volume factor ( |l5[ ) into (^), the cumulative distribution function is analytically 
evaluated as follows: 
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with the function V(a,x) being the incomplete Gamma function, F(a,x) = J°° dtt a ~ 1 e~ t . Here, the function i/r(<5) is defined 
by 

te(5)sa Ud Ia-J = a U) ' (19) 

where 5» = <$o(Rth, M*) and denotes the characteristic smoothed density of the halo with mass M*, which play a key role to 
reveal the sensitivity of the non-Gaussian tails. 

Equation ( |r?| ) indicates that a shape of the PDF is determined by the combination of two contributions: One is the halo 
profile. The scaling S~ 3 ^ a in equation ([171) reflects the slope of the halo profile. The flatter the halo profile is, the more rapidly 
the volume factor varies, thus the steeper the PDF becomes. Another contribution is the shape of the mass function, since 
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Figure 1. F(8, a, n) defined by eq.(^) as a function of 8/8*. Dependence of the halo density profiles {left) and initial spectra(rig/i4). 



the factor M/pA v i r gives the characteristic volume of the halo with the mass M, and the lower limit of the integration in 
equation (^) is, roughly speaking, determined by the halo mass whose characteristic smoothed density (defined by eq. jl3|]) 
just coincides with 8. Therefore, the function G(8; a, n) gives the volume fraction of halos (i.e., a sum of characteristic volume 
of halos) whose mass is larger than such lower limit, and directly reflects the shape of the halo mass function. Note that the 
dependence of the spectral index n on the PDF enters only through the function G(8; a, n). 
Now from (|l7|), the one-point PDF P(8;Rth) is analytically expressed as follows: 
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F(5;a,n) = 
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The expression @ with © and @ is a main analytical result in this paper. Qualitatively, it reveals the strong dependence 
of the halo density profiles and the sensitivity to the initial spectra. 

Figure [l] shows the behavior of the function F(S;a,n) for various values of n (left panel) and for a (right panel). Note 
again that the dependence of the power index n on the PDF enters only through the function F(S; a, n). It is clear from both 
panels that the slope of the function F becomes sensitive to both n and a in the range 8/8* > 1, corresponding to the mass 
range M > M» . The sensitivity to the power index n simply results from the steep exponential tails of the halo mass function, 
which becomes very sensitive to the spectral index in such mass range. And the sensitivity to the slope of the halo profile 
a is also explained from the exponential tails of mass function, that is, the variation in the relation between characteristic 
smoothed density So (eq. [Q) and halo mass M due to a small change of a is significantly enhanced through such exponential 
dependences. On the other hand, the slope of the function F, in turn, becomes insensitive to both n and a at the region 
8/8.* < 1. This is because in such range (corresponding to the mass range of M < M») the mass function only weakly depends 
on the spectral index n and scales as n(M) oc m'" -9 '/ 6 . Therefore, within the range 8/8* < 1, the PDF is almost solely 
determined by the contribution from the halo profile and varies as oc <5~' 3/ ' a+1 - ) , while at the region 8/8* > 1, the function 
F varies much faster than J _ ( 3 / a+1 ) ] thus the shape of PDF is determined by the contribution from the exponential tails of 
mass function. 

Figure |^ shows the approximate expression of one-point PDF ( pp| ) for various slope of the halo profiles (left) and for various 
initial spectra (right). For definiteness, we assume the Einstein-de Sitter universe (density parameter Qq — 1, cosmological 
constant Ao = 0) and the non-linear mass M* in each model is determined by introducing the artificial length scale Lbox an d 
normalizing the rms fluctuation a M to unity at Rm = O.lLbox- The smoothing radii are specifically chosen as both the large 
8* (~ 30 — 90) (Rth = 0.02Lbox) and the small 8 (~ 0.3 — 1.5) (Rth = 0.15Lbox)- Table Q summarizes the values of 8* for the 
models shown in Figure ^. 

In left panel, the spectral index are fixed to n = — 1 but the slope of halo profile a varies from 1 to 2.5, while in right 
panel, a is fixed to 2 but n varies from —2 to +1. It is clear from the cases of i? t h = 0.02ib O x in both panels that at the 
region 8 < 8*, the slope of PDF hardly depends on n but sensitively depends on a. On the other hand, at 8 > 8* (see the 
cases of Rth = 0.02ib O x in the right panel), the shape of PDF strongly depends on n. It follows from these results that if one 
takes smaller smoothing scale than the virial radius of the halo mass M* (thus 8* > 100), the slope of the PDF almost solely 
determined by the slope of the halo profile at the intermediate range of 5 < 8 < 50, while the shape of PDF strongly depends 
both on a and n if we conversely take a larger smoothing scale (so that 8* < 10). Figure ^| shows this point clearly. In the 
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Fig ure 2. Dependence of the halo density profiles (left) and initial spectra(ri<?/ii) on the approximate expression of the one-point PDF 
(|20|) in scale-free models. Here, the normalization of rms mass fluctuation <tm is determined by introducing the auxiliary length scale 
Lbox an d setting Gyi = 1 at Rm = O.lLbox m each panel. While the initial spectrum and the top-hat smoothing radius are respectively 
chosen as n = — 1 and R t ^ = 0.02Lbox m left panel, the right panel fixes the slope of the power-law profile to a = 2.0. 



Table X. The characteristic smoothed density 8* for models presented in Figure 
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figure, the PDF in the case of n = — 1 and a — 2 is shown restricting the contribution of halo mass to some small mass ranges. 
It turns out that at 8 > 8,, the major contribution to the PDF at a certain value of 8 comes from halos whose characteristic 
mass scale So roughly coincides with 8, thus the shape of PDF reflects shape of the halo mass function. At 8 < 8*, however, 
halos with So = 8 no longer lead to a major contribution, but halos whose characteristic density So is larger than 8 instead 
contribute to the PDF. In this case, PDF reflects the profile of the halos. 



4 COMPARISON WITH N-BODY SIMULATIONS 

We are now in a position to compare the one-point PDFs from the dark halo approach with those obtained from N-body 
simulations. To investigate this issue, the analytic results in the previous section turn out to become inadequate because of the 
non power-law nature of the realistic halo profiles. Hence, all the model predictions presented below are based on the NFW 
profile (Q). In computing the analytic PDF ([il^), the top-hat smoothed halo p 3m ooth is numerically evaluated in the Fourier 
space, in which we use the following fitting function y(k; M) for the Fourier transform of the NFW profile p N Fw/p (Ma & Fry 
2000): 

v(k . M ) - ^8 r 3 Me + l/g)-ln[hi(e + l/g)] m) 
with q = kr B . 

For the present purpose, we specifically treat the N-body simulation data for the scale-free models with initial power 
spectra P(k) oc k n (n = 1, 0, — 1, — 2)(Jing 1998) and the CDM model with cosmological constant(Lambda CDM, hereafter 
LCDM; Jing & Suto 1998). All the models employ N = 256 3 dark matter particles in a periodic comoving cube, L 3 ox , where 
the box size of the LCDM model is specifically chosen as Lbox = 100/i _1 Mpc. The gravitational force calculation is based 
on the P 3 M algorithm. While the scale-free models assume an Einstein-de Sitter universe, cosmological parameters of the 
LCDM model are chosen as (fio, Ao, h, as) = (0.3,0.7,0.7, 1.0), where the normalization as means the linear rms fluctuation 
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Figure 3. Contribution to the tails of PDF from different halo mass range for the model with n = —1 and a = 2. The thick solid lines 
are the exact calculation of equation ( Jlo[ ) with (^), while the thin lines represent the results restricting the range of integration. Here, 
the smoothing radius is specifically chosen as R t h = 0.02[Li DOX ], but the the horizontal scaling (5/(5* is not affected by any specific choice 
of flth- 



Table 2. The estimated values of the rms fluctuation amplitude cr(i?th) m Figure [j from the simulation data (labeled cr a ; m , the non-linear 
fitting formula by Peacock & Dodds (1996) (labeled by opd) and (5*. 
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at i? t h = 8/i _1 Mpc. As for the scale- free models, the normalization of the density fluctuation is determined by setting the 
linear rms fluctuation to unity at Rth = O.libox- 

Figure [| shows the N-body results of the one-point PDF for the top-hat smoothed density field in various smoothing radii 
(crosses). The left panel shows the results in LCDM model, while the right panel plots the PDFs in the scale-free model with 
the initial spectrum n = — 1. The error-bars in the left panel indicate the 1-cr errors among three different realizations. The 
smoothing radii in left(right) panel are chosen as R t h = 2, 6, 18/i _1 Mpc (Rth = 0.02, 0.05, 0.15Z/b ox ),which are typically scaled 
from non-linear to weakly non-linear regime. In Table ^, computed values of the rms amplitude a are presented (labeled by 
o"si m ). Obviously, the degree of non-linearity becomes significant as decreasing Rth and the resultant PDFs show much longer 
non-Gaussian tails. 

In Figure ^, the one-point PDFs from the dark halo approach are shown in solid lines. Within the validity range S ^ 1, 
in both panels, the dark halo approach reproduces the non-Gaussian tails of simulated PDF reasonably well. The predictions 
at large smoothing radius R t h = 18/i~ 1 Mpc(0.15I/box) somewhat over-predict the simulations, however, as have been pointed 
out by Kayo et al.(2001), this is most likely due to the finite volume effect of the N-body simulations. 

For comparison, the dashed lines in Figure ^ show the empirical log-normal distribution. The one-point PDF of the 
log-normal model is analytically expressed as 

P (X\ - 1 r-:n f + S) + aj N /2] 2 \ 1 

Tln(o) - exp < f 7~rr> i 23 ) 

where the single parameter ctln is related to the rms amplitude of density fluctuation at a given smoothing radius R t h, o~(Rth)' 

at N =log[l + ff 2 (ilt h )]. (24) 

For a given set of cosmological parameters, the quantity a (Rth) is computed using a fitting formula of non-linear power 
spectrum by Peacock & Dodds (1996) or even by the dark halo approach (Seljak 2000; Ma & Fry 2000; Scoccimarro et al. 
2001; Hamana, Yoshida & Suto 2001). The resultant values are summarized in Table (labeled by <tpd) and adopting those, 
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Figure 4. Comparison of one-point PDF between analytic models (solid, dashed) and N-body simulations (crosses) in LCDM model(£e/t) 
and scale-free model with n = —l(right). In eac h p anel, the results from dark halo approach are shown in solid lines, while the dashed 
lines indicate the log-normal model prediction In plotting the log- normal PDF, the rms fluctuation of density field a is evaluated 

using the fitting formula by Peacock & Dodds (1996). The error-bars in left panel denote the 1-cr errors among three different realization 
of simulation data. Also, in the right panel, the asymptotic slopes of the dependence £ — ("/S+l) with a = 1, 2, and 3 are shown. 




1 10 1 10 10 100 10 100 
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Figure 5. Comparison of non-Gaussian tails of PDF between dark halo approach and N-body simulations in cases with scale- free initial 
spectra. The left(right) panel depicts the results at the smoothing radius R t ^ = 0.02Lbox(^th = 0.15Lb ox ). F° r comparison, we also 
plot the log-normal PDFs in the right-side of the right panel (thin lines), where rms fluctuation of log-normal PDF, <r is calculated from 
simulations. Also, the asymptotic slope of the dependence <5 — ( c "/ 3 + 1 ) with a = 1, 2, and 3 are shown. 



the log-normal predictions ( |23| ) are plotted. In this sense, both empirical and analytic models in Figure ^ provide a fairly 
self-consistent prediction without invoking any information of simulation data. 

As mentioned in Kayo et al. (2001), the rms amplitude upd takes a slightly larger value than that of the simulation, <j s im, 
which should be ascribed to the finite box-size of the simulations. Apart from this small discrepancy, the log-normal PDFs 
approximate the simulated PDF fairly accurately in both low-density and high-density regions. The most striking evidence 
is that the tails of analytic PDF almost resembles those of the log-normal PDF, which is also consistent with the prediction 
from perturbation theory (Bernardeau 1994b; Bernardeau & Kofman 1995). According to the Table ^, the theoretical estimate 
of characteristic smoothed density 8, roughly varies from <5* ~ 0.03 to 40 as decreasing the smoothing radius. Thus, the 
tails of PDF at smoothing radius i? t h = 2/t -1 Mpc or 0.02Lb ox basically reflect the shape of the halo profiles and they scale 
as <5~( a / 3+1 ). Comparing the shape of PDFs with the slope of the dependence <5 _ ( a / 3+1 ) superimposed in the right panel of 
Figure ^, the non-Gaussian tails at smaller smoothing radii can be regarded as a consequence of outer slope of the NFW 
profile, a — 3. 

Figure [| shows weak and strong sensitivity of the non-Gaussian tails to the initial spectrum in scale-free models. Specifi- 
cally choosing the smoothing radius as Rth = 0.15Lb O x('e/¥) and R t h = 0.02Lbox(right), the results from the dark halo approach 
are plotted in right-side of each panel, while the left-side shows the simulation results. For comparison, the log-normal PDFs 
are also plotted in right-side of the right p&ne\(thin-lines) , where rms amplitude a is calculated from simulations. In cases with 
large smoothing radius Rth = 0.15Lb O x, the simulated PDFs rapidly fall off and they eventually vanish, while the analytic 
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model of PDF shows a slightly longer non-Gaussian tail. This discrepancy is basically ascribed to the finite volume size of 
the simulations (Kayo et al. 2001), however, the tendency seen in both panels is qualitatively similar. The analytic model 
reproduces the systematic behavior for the variation of spectral index n (see also the right-panel of Fig.^|). On the other hand, 
for smaller smoothing radius with R t ^ — 0.02Lb ox , except for the noisy data of scale-free model with n = +1, both left- and 
right-side of the panel show not only qualitatively but also quantitatively similar dependence on the initial power spectra. 
Again, the asymptotic behavior of the PDF is quite similar to the results with steeper power-law profile a ~ 3 (compare 
with the slope 5~( a / 3+1 ) superimposed in the figure). As a consequence, the non-Gaussian tails of PDF seem to resemble the 
log-normal PDF insensitive to the initial spectrum, as has been reported by Kayo et al. (2001). At a closer look at tails of 
PDF, however, the simulated PDFs still exhibit the initial spectrum dependence and they slightly deviate from log-normal 
PDFs (thin-lines). In contrast, the analytic PDF from the dark halo approach successfully reproduces this weak dependence. 
Note that the theoretical estimate of the characteristic smoothed density yields 5* ~ 10-70 at Rth = 0.02ibox, which suggests 
that the the simulated PDF will show a strong sensitivity to the spectral index at the high density tail 8/8, 3> 1, and thereby 
a large discrepancy between log-normal prediction and simulation is expected. Interestingly, there is some indication at the 
region 8 £ 200 in both simulation and halo approach, where the systematic variation for n apparently ceases. 

Therefore, we conclude that the analytic model based on the dark halo approach reasonably predicts the general tendency 
in both weakly and strongly non-linear regime and it quantitatively provides a useful approximation to the non-Gaussian tails 
of simulated PDF in the strongly non-linear regime. 



5 DISCUSSION & CONCLUSION 

In this paper, utilizing the dark halo approach, we have presented a simple analytic model characterizing the non-Gaussian 
tails of one-point PDF and investigated the sensitivity to both the initial power spectrum and the halo density profiles in the 
non-linear regime. Assuming the power-law profile of dark halos and the scale-free initial spectrum, the approximate expression 
for the tails of PDF is analytically derived (eq.po| with jl9j [^l|). Using this result, we explained the sensitivity to the initial 
spectrum seen in the simulation and the significant dependences on the halo profiles. Further, quantitative comparison with 
N-body simulation was made, which leads to the conclusion that the analytic model reasonably reproduces the simulated 
PDF qualitatively and quantitatively in a strongly non-linear regime. Then, the tails of simulated PDF turn out to resemble 
the log-normal distribution ( p3| ) at the region 8 ^ 100 in a strongly non-linear regime. 

In our model, the key quantity is the characteristic smoothed density 8, = p smoo th(0; R, M* )/~p, which discriminates the 
weak/strong sensitivity to the halo profile and initial power spectrum (see eq.(lS[]). At 8/8, 1, the tails of PDF basically 
follow the halo profile and varies as 5~^ 3 ' a+1 ', while the PDF at 5/8* ^ 1 exhibits a strong sensitivity to the initial spectrum, 
which can be ascribed to the steep exponential tails of halo mass function. In CDM models at present, the characteristic 
smoothed density typically takes the value 5, 3> 1 for a small smoothing scale Rth ~ lMpc//i and conversely i, < 1 for a 
large smoothing scale R t h > 10Mpc/h. Therefore, the log-normal behavior seen in the non-linear regime can be regarded as a 
consequence of the universality of the halo density profile. Further considering the fact that the slope of the simulated PDF is 
similar to that of the analytic results with a ~ 3, the weak sensitivity to the initial spectrum is especially due to the universal 
outer profile of p oc r -3 . 

The qualitative feature in the dark matter PDFs might be also the case for the PDF of galaxies, which has been known 
to be well approximated by the log-normal PDF (e.g., Hamilton 1985; Bouchet et al. 1993; Kofman et al. 1994) as the density 
PDF is. The number density of galaxies within the clusters is known to be fitted well by the King model (King 1962), which 
has a radial profile of p g oc [1 + (r/r c ) 2 ] -3 ^ 2 , where r c is a core radius of order 0.1/i -1 Mpc (Wu & Hammer 1993 and references 
therein). Thus, on suitably larger smoothing scale _R t h ~ r c , smoothed density profile becomes p g oc r -3 , at an outer part. 
Therefore, it is very likely that the similarity in the shape of the PDF tails between the dark matter and the galaxy is due to 
simply their similar distribution within halos. 

Further, the good agreement in the non-Gaussian tails of density PDF between the halo approach prediction and N-body 
data indicates that the halo model is not only used to qualitative study but also allows us to make a quantitative prediction. 
Although the dark matter PDF is hard to measure, it may be applied to related measurements such like the weak lensing 
convergence which is a line-of-sight projection of the matter density weighted by the lensing efficiency and is obtained from 
a coherent distortion in distant galaxy images (see Mellier 1999; Bartelmann & Schneider 2000 for reviews). Taruya et al. 
(2002) showed that the lensing convergence PDF is well approximated by the log-normal model and pointed out that the 
lensing convergence field can be regarded as a simple projection of the density field. This implies that one can safely apply the 
halo approach to the lensing convergence PDF. Actually, Kruse & Schneider (2000) employed the halo approach to predict 
the non-Gaussian tail in PDF of lensing convergence arguing it may be used to discriminate the cosmological model. Their 
prediction was tested against numerical experiments and a good agreement was found (Reblinsky et al. 1999). Since the 
approach adopted by Kruse & Schneider (2000) is essentially the same as described in this paper, one naturally expects that 
the tails of lensing convergence PDF basically follows the outer slope of the halo profiles at smaller angular scale 9 ~ l'-5', 
and the log-normal behavior found by Taruya et al. (2002) can be explained by similar reason to the mass density PDF. 

Therefore, as far as the CDM model with a Gaussian initial condition is concerned, theoretical prediction based on the 
dark halo approach consistently explains the log-normal nature of observable clustering statistics as a common feature in the 
hierarchical clustering universe. This conversely indicates that the statistics related to the non-Gaussian tails of PDF provide 
a powerful cosmological probe to the nature of dark matter clustering, which enables us to constrain the number density and 
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the density profile of dark halos. Probably, one promising application is the weak lensing statistics, as has been demonstrated 
by Kruse & Schneider (2000). Recently, from the 2.1 deg 2 image taken with a wide field camera on the prime focus of Subaru 
Telescope, non-Gaussian signature of the lensing convergence was detected evaluating the peak distribution function (Miyazaki 
et al. 2002). The number count of high positive peaks, which are thought to come from the dark halos, was compared with 
the theoretical prediction using halo approach and showed a consistent result with NFW profiles under the LCDM cosmology. 
Note, however, that the analytic models of PDF from the dark halo approach currently limit the applicability to highly non- 
Gaussian tails in a strongly nonlinear regime. To obtain a more stringent constraint on the cosmological parameter and/or 
the number density of dark halos, we should further improve the applicability to the low-density region up to the observable 
scale of the non-linear regime. In the light of this, the perturbative construction of PDF by Bernardeau & Valageas (2000) 
or the non-perturbative modeling from the hierarchical ansatz (Munshi & Jain 2000; Valageas 2000), as well as the empirical 
log-normal model (Taruya et al.2002) could be useful and be regarded as a complementary approach. Combining the present 
work with these approaches would be possible to improve the prediction dramatically and will be discussed elsewhere. 
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APPENDIX A: SMOOTHED DENSITY PROFILE FOR POWER-LAW HALOS 



In this appendix, adopting the top-hat filter function we derive the exact expression for the smoothed halo profiles. 

Substituting the definition ([H]) into (j^), the expressions for smoothed density profile can be reduced to 



,(r;R th ,M) = 



Jo 



dV e{R th -\r'\) p(\r-r'\;M) 



ds s 



dx p(\f: 



r 1 + s 2 



2rsx ; M). 



(25) 



For the power-law profiles M, all the integrals in the expression (B5J) are analytically evaluated, which lead to 



oth (r;7? th ,M) = 



36 c pr° 
2i? 3 



2 (a - 2) 



ds s 



2 r- 



(r + s) 2 



Ppowcr(-R th ;M) 

3M(i--) 3 ~' 



(a — 2)r s 
+ (l+u)'^} + ^{(l-u) i 

-(u-l) 3 -} + 



(1 + u) 4 -} 

(1+u) 4 -*} 



(0 < r «S R t 
(Rth < r) 



(26) 



Here, the variable u denotes r/Rth- Notice that the above expressions are valid for the profiles with fl^2. For the power- law 
profile with a — 2 (singular isothermal sphere), the logarithmic dependence of the radius r appears, which can be separately 
evaluated as follows: 



oth(r;i?th,M) = 



3Sc prj 
2^4 



ds — 



r + s 



\p{R^,M) ji+ifu-i-u) logj^j}. 



(27) 



While the above results are exact, it is convenient and useful to derive the simple analytic expressions for the one-point PDF 
P(S\ Rth) (see Sec.^). From the asymptotic behavior of the smoothed halo profiles, we approximate the expressions ( |26[ ) and 
(p7|) as follows: 



(3=5) P P owc r (i? th ;M) 



oth(r;7?th,M) 



l/a 



1/a 



Rt] 



Rth 



(28) 



The left panel of Figure ^ compares the approximation ( pq ) with the exact results ( |26[ ) and (|27j). Apart from the intermediate 
scales r/Rth ~ 1, the approximation recovers the exact results reasonably well. On the other hand, deviation at the intermediate 
scale slightly increases as increasing the slope of profile a, which leads to the under-prediction of volume factor V(> 5; R t h,M). 
Hence, for steeper halo profile a ~ 3, the analytic PDF using (^i|) tends to under-predict the amplitude of PDF as shown in 
the right panel. 
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Figure 6. Left panel: comparison of smoothed power-law profiles between exact results (Eqs.|2fj] J27|; solid) and approximations (Eq:J28|; 
dashed). Right panel: the PDFs computed from the exact results without using the approximation ( p^ ) (thin lines) compared with the 
approximation results. 



To check the validity of the approximation ( J2C| ) , together with the analytic PDF in the left panel of Figure EL the PDF 
computed from the exact results without using approximation ( |l2| ) is plotted and the results are presented in the right panel 
of Figure ^| (thin lines). As is expected from the left panel, the approximate expression ( po[ ) for steeper halo profiles with 
q = 2.5 tends to lower its amplitude, but the qualitative behavior still remains correct. 



APPENDIX B: EFFECTS ON UNCERTAINTY OF THE MODEL PREDICTION 

The analytic model of one-point PDF presented in section ^] suffers from uncertainties in the model assumptions (mass function 
of dark halos and halo density profile). Among these, the concentration parameter of the NFW halo profiles, c(M) might 
be crucial in comparing the analytic prediction with simulations (Sec||), since the sensitivity to c(M) has been reported in 
previous study using dark halo approach (e.g., Seljak 2000). Although we adopt the fitting form of the concentration parameter 
by Bullock et al. (20001) with a particular set of the parameters (c„, /3) (see eq. Q below), we examine the influence of varying 
the parameters (c»,/3). 

Figure |?] shows the dependence of the non-Gaussian tails on the variation in concentration parameter assuming the LCDM 
universe. Upper(lower) panels represent the results fixing /3(c*). Although the influence of the variation (c,,/3) turns out to be 
small, the systematic behavior in the tails of PDF seems to vary in a opposite sense: decreasing the parameter j3 or increasing 
c„ suppresses the non-Gaussian tails. 

The small influence on the tails of PDF can be easily deduced from the Fourier transform of the NFW profile (^2|), in 
which the concentration parameter appears only through the factor 5 c rl, leading to the logarithmic dependence of c(M). On 
the other hand, apparently opposite behavior for the variation of (c»,/3) is explained as follows. First note that decrease of (3 
or increase of c* yields more concentrated density field for halos more massive than the nonlinear mass A/* . This implies that 
while the number of halos in which the central density exceeding a specific value of S increases, the volume of each density 
profile exceeding S conversely decreases in the massive halos. Since the contribution of massive halo is relatively dominant in 
the tails of PDF, the latter effect finally wins and accordingly the amplitude of one-point PDF decreases. 

Except for the details, the systematic uncertainty of the concentration parameter does not affect the final conclusions in 
section ^ and the dispersion of 0.2 in lnc (Jing 2000; Bullock et al. 2001) can be safely neglected in our calculation. Note, 
however, that in the case of the lensing convergence PDF, the dependence of (c, ,(3) becomes slightly larger than that in the 
mass density PDF and the systematic behavior of non-Gaussian tails appears in a opposite sense to P(S). These differences 
simply result from the projection of the three-dimensional density field to the two-dimensional field (e.g., Takada & Jain 
2002). 
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